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RELIABLE IDENTIFICATION OF COMPTON-THICK QUASARS AT Z ; 
SPITZER MID-INFRARED SPECTROSCOPY OF HDF-OMD49 
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ABSTRACT 

Many models that seek to explain the origin of the unresolved X-ray background predict that Compton-thick 
Active Galactic Nuclei (AGNs) are ubiquitious at high redshift. However, few distant Compton-thick AGNs 
have been reliably identified to date. Here we present Spitzer-TRS spectroscopy and 3.6-70 photometry of a 
Z = 2.211 optically identified AGN (HDF-oMD49) that is formally undetected in the 2 Ms Chandra Deep Field- 
North (CDF-N) survey. The Spitzer-TRS spectrum and spectral energy distribution of this object is AGN domi- 
nated, and a comparison of the energetics at X-ray wavelengths to those derived from mid-infrared (mid-IR) and 
optical spectroscopy shows that the AGN is intrinsically luminous (Z-2-10 keV ~ 3 x 10 44 erg s _1 ) but heavily ab- 
sorbed by Compton-thick material (7Vh 3> 10 24 cm -2 ); i.e., this object is a Compton-thick quasar. Adopting the 
same approach that we applied to HDF-oMD49, we found a further six objects at z w 2-2.5 in the literature that 
are also X-ray weak/undetected but have evidence for AGN activity from optical and/or mid-IR spectroscopy, 
and show that all of these sources are likely to be Compton-thick quasars with L2-i0keV > 10 44 erg s -1 . On 
the basis of the definition of Daddi et al. (2007), these Compton-thick quasars would be classified as mid-IR 
excess galaxies, and our study provides the first spectroscopic confirmation of Compton-thick AGN activity 
in a subsample of these z ~ 2 mid-IR bright galaxies. Using the four objects that lie in the CDF-N field, 
we estimate the space density of reliably identified Compton-thick quasars [<t> « (0.7-2.5) xl0~ 5 Mpc~ 3 for 
^2-iOkeV > 10 44 erg s objects at z ~ 2-2.5] and show that Compton-thick accretion was probaby as ubiqui- 
tious as unobscured accretion in the distant Universe. 

Subject headings: galaxies: active — galaxies: high-redshift — infrared: galaxies — X-rays: galaxies — 
ultraviolet: galaxies 



1. INTRODUCTION 

There is a growing need for a complete census of Active 
Galactic Nuclei (AGNs). The seminal discovery that every 
massive galaxy in the local Universe harbors a super-massive 
black hole (Mbh ~ 10 6 M ) implies that all massive galaxies 
must have hosted AGN activity at some time during the past 
w 13 Gyrs (e.g., Rees 1984; Kormendy & Richstone 1995). 
To trace accurately how and when these black holes grew re- 
quires a detailed census of AGN activity that will provide, 
amongst other things, estimates of the efficiency and duty cy- 
cle of black-hole growth and any dependencies of nuclear 
obscuration on AGN luminosity, redshift, and host-galaxy 
type (e.g., Marconi et al. 2004; La Franca et al. 2005; Shankar 
et al. 2007). 

The exceptional sensitivity of the Chandra Deep Field ob- 
servations (e.g., Brandt et al. 2001; Giacconi et al. 2002; 
Alexander et al. 2003a) has helped to unveil a > 10 times 
larger population of AGNs than found at most other wave- 
lengths (« 7200 deg~ 2 ; e.g., Bauer et al. 2004). Optical 
spectroscopic follow-up observations have shown that these 
AGNs are detected out to z ~ 5, and detailed X-ray spectral 
analyses have revealed that the majority of the sources are 
obscured by gas and dust (e.g., Barger et al. 2003; Szokoly 
et al. 2004; Tozzi et al. 2006; see Brandt & Hasinger 2005 



for a review). However, although clearly effective at identi- 
fying even heavily obscured AGNs out to high redshift, there 
is compelling evidence that a large fraction of the AGN pop- 
ulation remains undetected at the < 10 keV observed-frame 
energies probed by these surveys: (1) about half of the X-ray 
background (XRB) is unresolved at > 6 keV (Worsley et al. 
2005), (2) the observed obscured:unobscured AGN ratio is 
lower than that found for comparably luminous AGNs in the 
local Universe (e.g., Treister & Urry 2005), (3) the 5-10 keV 
blank-field number counts are steeply rising at the faintest X- 
ray fluxes (Rosati et al. 2002), and (4) few Compton-thick 
AGNs (TVh ^ 1.5 x 10 24 cm~ 2 ) have been identified, even 
though they comprise ks 50% of the AGN population in the lo- 
cal Universe (e.g., Risaliti et al. 1999; Guainazzi et al. 2005; 
Tozzi et al. 2006). 

Many of the X-ray undetected AGNs are expected to be 
intrinsically luminous sources that are heavily obscured by 
Compton-thick material (i.e., 7Vh > 1.5 x 10 24 cm~ 2 ; see Co- 
mastri 2004 for a review). The most robust identification of 
a Compton-thick AGN is made from high signal-to-noise ra- 
tio (S/N) X-ray spectroscopy, where the detection of a high 
equivalent-width Fe K emission line (Wx > 1 keV at rest- 
frame energies 6.4-6.9 keV) and a steeply rising reflection 
component at > 10 keV reveals that little or no X-ray emission 
is being seen directly, implying that the central source is very 
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heavily absorbed (e.g., George & Fabian 1991; Matt, Brandt, 
& Fabian 1996; Maiolino et al. 1998; Matt et al. 2000). How- 
ever, since the extreme obscuration toward the nucleus of a 
Compton-thick AGN renders the observed < 10 keV emis- 
sion ss 30-1000 times fainter than the intrinsic (i.e., unab- 
sorbed) emission, it is often challenging to identify robustly 
these sources on the basis of their X-ray data alone; the range 
of absorption correction factors is based on the observed to 
intrinsic X-ray luminosity ratio for the AGNs in Table 8.1 of 
Comastri (2004), with the observed X-ray fluxes from Bas- 
sani et al. (1999). For example, despite Compton-thick AGNs 
comprising half of the AGN population in the local Universe, 
only « 50 sources have been reliably identified from X-ray 
data to date (e.g., Comastri 2004). Fortunately, although not 
as conclusive as high S/N X-ray spectroscopy, the presence 
of Compton-thick absorption can also be identified in X-ray 
weak AGNs when a reliable measurement of the intrinsic 
power of the AGN is available at other wavelengths; we stress 
here the necessity for X-ray constraints in order to identify a 
Compton-thick AGN since it is only in the X-ray band that an 
absorbing column density can be measured or inferred. 

Two of the most reliable measurements of the intrinsic 
power of an obscured AGN are the luminosity of the mid- 
infrared (mid-IR; rest-frame > 3 pm; e.g., Efstathiou & 
Rowan-Robinson 1995; Granato et al. 1997; Lutz et al. 
2004) continuum and high-excitation emission lines (i.e., 
[OIII]A,5007; e.g., Allen 1973; Kwan & Krolik 1981; Bassani 
et al. 1999), both of which are believed to be directly illumi- 
nated by the central source. These emission regions provide 
a good proxy for the intrinsic AGN luminosity, even in the 
presence of extreme Compton-thick absorption, since they are 
more extended than the X-ray absorbing material (i.e., larger 
than the broad-line region; e.g., Lamer et al. 2003; Risaliti 
et al. 2007). As the luminosities of both the mid-IR contin- 
uum and the high-excitation emission-line region are depen- 
dent on many factors, including the location and geometry of 
the region with respect to the central source, more robust con- 
straints are placed if both measurements are available. 

Comprehensive evidence for a large X-ray undetected AGN 
population has been found from a variety of X-ray based anal- 
yses using sources detected in deep Spitzer surveys at mid-IR 
wavelengths, suggesting that Compton-thick AGNs may be 
ubiquitious in the distant Universe (e.g., Donley et al. 2005, 
2007; Alonso-Herrero et al. 2006; Polletta et al. 2006; 
Daddi et al. 2007a; Steffen et al. 2007; Martmez-Sansigre 
et al. 2007; Fiore et al. 2008). However, all of these stud- 
ies have relied on photometric observations for the identifi- 
cation of candidate Compton-thick AGNs, leading to signifi- 
cant uncertainties in measurements of the intrinsic AGN lumi- 
nosity, the absorbing column density, the degree of contami- 
nation from star-formation activity, and by implication, the 
space density of distant Compton-thick AGNs. In this paper 
we present Spitzer-IRS spectroscopy and 3.6-70 pm obser- 
vations of an optically identified narrow emission-line AGN 
at z = 2.211 (HDF-oMD49; Steidel et al. 2002) that is un- 
detected in the published catalogs of the deepest X-ray ob- 
servation currently available [the 2 Ms Chandra Deep Field- 
North (CDF-N); Alexander et al. 2003a]. From a compar- 
ison of the AGN energetics at X-ray wavelengths to those 
determined from optical and mid-IR spectroscopy, we ro- 
bustly show that HDF-oMD49 hosts an intrinsically luminous 
(^2-iOkeV ~ 3 x 10 44 erg s -1 ) Compton-thick AGN; i.e., this 



source is a Compton-thick quasar. We then use these X- 
ray-optical-mid-IR diagnostics to identify a further six dis- 
tant Compton-thick quasars in the literature and place limits 
on the space density of Compton-thick quasars at z « 2. We 
adopted H Q = 71 km s" 1 Mpc -1 , Cl M = 0.27, and £2 A = 0.73 
throughout. 

2. OBSERVATIONS 

HDF-oMD49 (optical position a 2 ooo = 12 h 37 m 04P34, 
5 2 ooo = +62°14'46"2; Steidel et al. 2002) was originally 
identified in the Lyman-Break Galaxy (LBG) searches of Stei- 
del etal.(2002; 2003). With a redshift of z = 2.211, HDF- 
oMD49 lies at the low-redshift end of the LBG population 
(median redshift z ~ 3). HDF-oMD49 is also identified as a 
BX/BM galaxy, based on the criteria of Steidel et al. (2004), 
and it is dubbed BM1156 in the BX/BM study of Reddy 
etal. (2006). The apparent contradiction in HDF-oMD49 be- 
ing selected as both an LBG and BX/BM is due to small dif- 
ferences in the photometric data used in these studies. 

2.1. Spitzer Data 

HDF-oMD49 is detected by IRAC at 3.6-8.0 ^m and MIPS 
at 24 pm in the Spitzer observations obtained as part of the 
Great Observatories Origins Deep Survey (GOODS) legacy 
project (P.I. M. Dickinson; R. Chary et al., in prep.); the 
3.6 pm, 4.5 pm, 5.8 pm, 8.0 pm, and 24 pm fluxes are 6.1 pJy, 
7.0 ply, 14.8 ply, 40.1 ply, and 380 ply, respectively. Al- 
though HDF-oMD49 is bright at 24 /jm, it is undetected in 
the ultra-deep 70 pm observations of Frayer et al. (2006; 
3 O upper limit of fiofim < 2.0 mJy), indicating that it has a 
"warm" mid-IR Spectral Energy Distribution (SED; e.g., Pa- 
povich et al. 2007). HDF-oMD49 is also detected at 16 pm 
with /i6,um = 222 ply using the "peak-up" imaging capability 
of Spitzer-IRS (H. Teplitz et al., in prep.). 

The Spitzer-IRS spectroscopy of HDF-oMD49 was ob- 
tained in spectral staring mode as part of program G02- 
20456 (P.I. R. Chary) using the Long-Low modules (LL; 
R w 64-128; Houck et al. 2004). HDF-oMD49 was observed 
with LL1 (20.5-37.9 pm) for 3 hrs, and additional observa- 
tions with LL1 (0.93 hrs) and LL2 (1.33 hrs; 14.3-21.2 pm) 
were obtained when HDF-oMD49 fell serendipitously into 
the slit of another target from the same program. We an- 
alyzed the 2D basic calibrated data files from the S 14.0.0 
Spitzer-IRS pipeline data; data reduction includes cleaning 
rogue pixels, fitting and removing the latent charge build-up, 
subtracting the sky, and averaging the 2D files together. Spec- 
tral extraction was performed using a 2 pixel wide window 
in SPICE and the data were calibrated using the same extrac- 
tion window on a standard-star spectrum. 2 The average RMS 
of the LL1 and LL2 spectra are 97 pjy and 120 ply, respec- 
tively. The Spitzer-IRS fluxes measured through the 16 pm 
and 24 pm filter spectral response curves are 180 ± 120 ply 
and 370 ±100 p}y, respectively, consistent with the observed 
fluxes in these bands. 3 For more details on the observations 

1 Here we define a quasar as an AGN with L2-l0kcV > >0 44 er g s _1 . This 
definition is consistent with that derived from the classical optical quasar 
threshold of Mb < —23 (e.g., Schmidt & Green 1983), assuming the (Xox re- 
lationship of Steffen et al. (2006) and the composite quasar spectrum of Van- 
den Berk etal. (2001); this definition is also consistent with Lx.* derived for 
X-ray unobscured AGNs (see Table 5 of Hasinger et al. 2005). 

2 See http://ssc.spitzer.caltech.edu/postbcd/spice.html for details of 
SPICE. 

3 The spectral response curves for the individual Spitzer instruments can 
be obtained from http://ssc.spitzer.caltech.edu/obs/ 




FIG. 1 . — Spitzer-IRS spectrum of HDF-oMD49 (solid line) for the LL1 (thick line) and LL2 (thin line) orders; the 1 o error arrays are plotted as a dotted line. 
Comparison spectra (a: the starburst galaxy M82 and the Compton-thick AGN NGC 1068; b: the composite spectrum of X-ray obscured Compton-thin quasars 
from Sturm et al. 2006) are shifted to z = 2.21 1 and overlaid on the data of HDF-oMD49. The mid-IR spectrum of HDF-oMD49 is consistent with that of an 
AGN; see §3.1 for quantitative constraints. 



and data reduction see Pope et al. (2008). 

2.2. Chandra Data 

HDF-oMD49 is undetected in the published 2 Ms Chandra 
catalogs of Alexander et al. (2003a). However, these cata- 
logs were produced using a conservative source-detection al- 
gorithm to minimize the number of spurious detections, and 
it is possible to identify robustly fainter sources when search- 
ing for X-ray emission associated with a known source (e.g., 
see §3.4.2 of Alexander et al. 2003a and §5 of Alexander 
et al. 2001). Indeed, adopting a WAVDETECT (Freeman et al. 
2002) false-positive probability threshold of 10~ 5 we detected 
significant X-ray emission (X-ray position a2ooo = 12 h 37 m 
04?32, 82000 = +62°14'46"4) within 0.25 arcsec of the opti- 
cal position of HDF-oMD49 in the 0.5-8 keV, 0.5-2 keV, and 
4-8 keV bands with 8.1 ±3.5, 6.7 ±2.8, and 4.9 ±2.6 counts, 
respectively (see also Laird et al. 2006). HDF-oMD49 re- 
mains undetected in all of the other four X-ray bands explored 
by Alexander et al. (2003a) and we calculated 3 O upper lim- 
its following §3.4.1 of Alexander et al. (2003a). Although the 
S/N of the X-ray data is low, the low background of Chan- 
dra allows sources with very low count rates to be reliably 
detected. For example, with a measured background in the 
detection aperture of HDF-oMD49 (« 2 pixel radius) of only 
w 1 counts in the 0.5-2 keV band, the probability of detecting 
w 7 counts by chance is w 10~ 5 ; furthermore, HDF-oMD49 
is detected in two independent bands, increasing the overall 
probability that the detection is real. 

The X-ray counts correspond to fluxes of 7.1 x 
10~ 17 erg s- 1 cm- 2 (0.5-8 keV), 2.6 x 10~ 17 erg s" 1 crrT 2 
(0.5-2 keV), and 1.5 x 10~ 16 erg s" 1 cirr 2 (4-8 keV), for 
r = 1.4 (i.e., the spectral slope of the XRB). The detection 
in the 4-8 keV band (rest-frame 12.8-25.7 keV) is partic- 



ularly notable since Chandra is significantly less sensitive 
at 4-8 keV than at < 4 keV energies, suggesting that HDF- 
oMD49 has a flat X-ray spectral slope; the non detection in 
the wider 2-8 keV band provides support for this hypothesis. 
The X-ray band ratio (the 2-8 keV/0.5-2 keV count-rate ra- 
tio) is < 1.9, which corresponds to F > 0.2 and is consistent 
with the X-ray spectral slope derived from the 4-8 keV/0.5- 
2 keV count-rate ratio (r « 0.4); the difference in the X-ray 
fluxes estimated using F = 0.4 instead of F — 1.4 is negligi- 
ble given the low S/N of the X-ray data ( ^ 15%). The cor- 
responding rest-frame 1.6-6.4 keV, 6.4-12.8 keV, and 12.8- 
25.7 keV luminosities (or 3 O upper limit) of HDF-oMD49 
are 10 42 erg s -1 , < 4 x 10 42 erg s -1 , and 6 x 10 42 erg s _1 , 
respectively. 

3. RESULTS 

3.1. The Spitzer SED 

In Fig. la we show the Spitzer-IRS spectrum of HDF- 
oMD49. The spectrum is noisy, particularly in the LL2 or- 
der. However, there is no significant emission from the dom- 
inant rest-frame 7.7 /jm Polycyclic Aromatic Hydrocarbon 
(PAH) feature found in star-forming galaxies, indicating that 
the mid-IR spectrum does not have a strong contribution from 
star-formation activity. By contrast, the Spitzer-IRS spectrum 
of HDF-oMD49 is similar to the nearby Compton-thick AGN 
NGC 1068 and even shows weak evidence for Si absorption 
at 9.7 /jm, a feature often seen in obscured AGNs (e.g., Shi 
et al. 2006; Hao et al. 2007; Spoon et al. 2007). The Spitzer- 
IRS spectrum of HDF-oMD49 is also qualitatively similar to 
the composite mid-IR spectrum of the X-ray obscured quasars 

A A I T? 

(Lx > 10 erg s ; Nh > 10 cm ) investigated by Sturm 
et al. (2006); see Fig. lb. 

We can quantify the contributions from AGN and star- 
formation activity to the mid-IR luminosity of HDF-oMD49 
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Table 1 . Robustly Identified z ~ 2 Compton-thick Quasars 







d L 


i-X.obs 




Lhya 


Lew 


i-Hell 




^■[OIII] 


Lx,6/im 


Lxfints 




Name 




(Mpc) 


(erg s~') 


(erg s~') 


(erg s~') 


(erg s~') 


(ergs -1 ) 


(erg s~') 


(erg s~') 


(erg s~') 


(ergs -1 ) 


Refs 


SMM J 123600+62 1047" 


2.002 


15750 


<42.4 


45.3* 












44.7 




1,2 


HDF-OMD49" 


2.211 


17800 


42.1 


45.2* 


42.5 


42.4 


41.9 


41.7 




44.6 


44.4 


3,4 


FSC 10214+4724 r 


2.285 


18550 


42.3 


44.9* 










42.6 


44.3 


44.4 


5,6,7 


SW J104406+583954 


2.430 


19990 


<43.4 


45.8* 


43.9 


43.5 


42.9 






45.2 


45.6 


8,9 


BX160" 


2.462 


20290 


<42.3 


44.9'' 


42.3 


41.8 








44.3 


44.1 


3,4 


BX1637" 


2.487 


20570 


<42.7 


44.8 <i 


42.3 










44.3 


44.1 


3,4 


SW J104409+585224 


2.540 


21100 


<43.4 


46A b 


43.8 


43.1 








45.9 


45.5 


8, 10 



Notes - The luminosities are given in logarithmic units of erg s . The X-ray luminosities are given in the rest-frame 2-10 keV band: Lx.obs is the observed 
X-ray luminosity, calculated from the observed-frame 0.5-2 keV luminosity assuming T = 1.4, Lx.6/jm is the X-ray luminosity implied from the rest-frame 6 jum 
AGN continumm, and Lx. lines is tne average X-ray luminosity implied from the emission-line luminosities; see §3.3 and §4.1 for more details. The references 
correspond to the X-ray, mid-IR, and optical emission-line data - 1: Alexander et al. (2005b); 2: Pope et al. (2008); 3: This paper; 4: Reddy et al. (2006); 5: 
Alexander et al. (2005a); 6: Teplitz et al. (2006); 7: Serjeant et al. (1998); 8: Polletta et al. (2006); 9: Weedman et al. (2006); 10: Polletta et al. (2008). 
"Object lies in the GOODS-N field. 

b Rest-frame 6 fjm luminosity is calculated from the best-fitting AGN component to the Spitzer-IRS data. 
c Data have been corrected for an assumed lensing boost of 50. 
d Rest-frame 6 fjm luminosity is calculated from the Spitzer SED. 



\ Rest-frame Wavelength (jum) 10 




Observed-frame Wavelength (/Am) 

FIG. 2. — Infrared data of HDF-oMD49, showing the Spitzer-IRS spec- 
troscopy (thick solid curve) and Spitzer photometry (solid circle) at 3.6- 
70 ;um. The error bar in the .v-axis direction indicates the width that corre- 
sponds to > 50% of the maximum transmission in each of the photometric 
bands; see Footnote 3. The data are compared to the SEDs of the Compton- 
thick AGN NGC 1068 (dashed curve) and the starburst galaxy M82 (long- 
dashed curve), which are shifted to z = 2.21 1 and normalized based on the 
results from fitting the mid-IR spectrum of HDF-oMD49 (see §3.1). In order 
to fit the data in the IRAC bands it was necessary to also add a black-body 
component of very hot 1000 K dust (dotted curve). The thin solid curve shows 
the combined model results. The data for the NGC 1068 and M82 SEDs were 
taken from Rigopoulou et al. (1999), Chary & Elbaz (2001), Galliano et al. 
(2003), and Magnelli et al. (2007). 



by fitting the Spitzer-IRS spectrum with AGN and star- 
formation templates. Adopting the approach of Pope et al. 
(2008), we fit the Spitzer-IRS spectrum with an M82 starburst 
component and include an AGN component using the nuclear 
spectrum of NGC 1068, performing a % 2 minimization to ob- 
tain the best-fitting normalization. On the basis of this ap- 



proach, we find that the mid-IR spectrum of HDF-oMD49 is 
dominated by AGN activity, with a limit of a 10% contribution 
from star formation in the Spitzer-IRS band; there is no clear 
signature for star formation and this limit effectively corre- 
sponds to the noise per pixel in the spectra. We get consistent 
results if we fit the mid-IR spectrum with an M82 starburst 
component and an absorbed power-law component (to repre- 
sent the AGN), leaving the slope of the power law and the ex- 
tinction as free parameters. The best-fitting results with this 
second approach give a spectral slope of a = 2 (consistent 
with that found for AGN-dominated sources; e.g., Alonso- 
Herrero et al. 2003, where f v « v~ a ) with 19.7pm « 0.65 and 
a « 5% limit on the contribution from star formation. 

In Fig. 2 we show the 3.6-70 pm SED of HDF-oMD49 and 
plot on the NGC 1068 and M82 templates, normalized to the 
values obtained from the mid-IR spectral fitting. The strong 
observed-frame 3.6-8 pm emission and the power-law-like 
featureless SED strengthens our conclusion from the mid-IR 
spectral fitting, showing that the IR SED of HDF-oMD49 is 
dominated by AGN activity. The strongly rising continuum at 
> 4.5 pm (rest-frame > 1.4 pm) also indicates that the con- 
tribution from starlight (which typically peaks at rest-frame 
1.6 pm) is weak at these wavelengths. The combination of 
the NGC 1068 and M82 templates gives a respectable fit to 
the overall SED of HDF-oMD49 even though the wavelength 
coverage of the SED data is three times broader than that used 
when fitting the mid-IR spectrum (rest-frame 1.1-21.8 pm 
versus rest-frame 4.4-10.6 pm). However, the combination 
of the NGC 1068 and M82 components lies below the ob- 
served emission in the IRAC bands (rest-frame < 3 pm). This 
excess emission can be fitted by adding a black-body com- 
ponent of very hot dust (« 1000 K), as often found in ob- 
scured AGNs (e.g., Alonso-Herrero et al. 2001). On the basis 
of the fitted NGC 1068 component, the rest-frame 6 lu- 
minosity from HDF-oMD49 is vL 6;Um w 1.3 x 10 45 erg s~ l 
(/6/im = 210 piy), indicating that the AGN is powerful (e.g., 
compare to the nearby AGNs explored by Lutz et al. 2004). 

Our analyses show that the IR SED of HDF-oMD49 is dom- 
inated by AGN activity and allow for just a small contribu- 
tion from star formation. Using the best fitting star-formation 
component normalizations and the relationship between PAH 
luminosity and star-formation rate (SFR) given by Pope et al. 
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(2008; see Eqn. 7), we predict SFR « 60-120 M yr" 1 , with 
the range representing the two different star-formation contri- 
butions we estimated from the mid-IR spectral fitting. The ra- 
dio detection of HDF-oMD49 with f\ 4ghz «32^Jy (G. Mor- 
rison et al., in prep.) implies a SFR of « 500 M @ yr _1 (calcu- 
lated following Bauer et al. 2002). This is clearly inconsistent 
with the mid-IR data, suggesting that either the radio emission 
is dominated by the AGN or the radio-SFR relationship is not 
applicable for high-redshift objects such as HDF-oMD49. 

3.2. The Chandra SED 

The contrast between the weak X-ray emission (see §2.2) 
and the bright AGN-dominated IR SED implies that the 
AGN in HDF-oMD49 is heavily absorbed at X-ray ener- 
gies. The signature of absorption should be evident in the 
X-ray spectrum of HDF-oMD49. However, before examin- 
ing the X-ray data it is necessary to consider the contribution 
to the X-ray emission from star formation. Assuming a SFR 
of 120 Mq yr (see §3.1) and the empirically determined 
Lx-SFR relationship found by Bauer et al. (2002), we pre- 
dict a rest-frame 0.5-8 keV flux of 2 x 10~ 17 erg s cm~ 2 , 
which corresponds to an observed-frame 0.5-2 keV flux of 
1.1 x 10~ 17 erg s~ ! cm~ 2 (Lx ~ 4 x 10 41 erg s _1 ), for a typi- 
cal X-ray spectral slope for star-forming galaxies of F = 1.8. 
Consistent X-ray luminosities are obtained from the empri- 
cially derived Ranalli et al. (2003) relationship, although the 
Persic et al. (2004) relationship would predict an w 4 times 
lower X-ray luminosity but only takes into account the contri- 
bution from high-mass X-ray binaries. These analyses show 
that the contribution from star formation could be significant 
to the observed-frame 0.5-2 keV flux (« 40%). However, 
the contribution will be negligible in the other X-ray band in 
which HDF-oMD49 is detected (observed-frame 4-8 keV). 

In Fig. 3 we show the X-ray data of HDF-oMD49 and 
compare them to the AGN model SEDs adopted by Gilli 
et al. (2007) to interpret the XRB. The X-ray data of HDF- 
oMD49 clearly indicate that the AGN is absorbed at X-ray 
energies; on the basis of the X-ray band ratio (see §2.2), the 
X-ray data are consistent with Na > 10 23 cm~ 2 (e.g., see 
Fig. 4 of Alexander et al. 2003b). However, with our low S/N 
X-ray data we cannot unambiguously distinguish between 
Compton-thick and Compton-thin absorption. Furthermore, 
the X-ray spectrum of a Compton-thick AGN is source spe- 
cific and depends upon the strength of the reflected and scat- 
tered X-ray components, in addition to the underlying power- 
law continuum and any additional absorption; for example, 
see Appendix A of Guainazzi et al. (2005) for information 
on the difficulty of identifying Compton-thick AGNs on the 
basis of X-ray colors alone. Therefore, without detailed X- 
ray spectral signatures that clearly indicate the presence of 
Compton-thick absorption (e.g., the identification of a promi- 
nent Fe K emission line or a strong reflection component), it is 
not possible to argue solely on the basis of low S/N X-ray data 
that an AGN is absorbed by Compton-thick material. In order 
to confirm whether such an X-ray weak source is Compton 
thick, it is necessary to estimate the intrinsic AGN luminos- 
ity. 

3.3. Quantifying the Intrinsic Properties of the AGN in 
HDF-oMD49 

A good proxy for the intrinsic luminosity of AGN activ- 
ity is the rest-frame mid-IR emission, which can provide 
an absorption-independent measure of the intrinsic luminos- 
ity even in the presence of Compton-thick absorption (e.g., 
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FIG. 3 .— X-ray data of HDF-oMD49, compared to the AGN model SEDs 
adopted by Gilli et al. (2007) to interpret the XRB (model Nn indicated) 
and normalized to the observed-frame 4—8 keV flux. The open circle shows 
the X-ray emission at observed-frame 0.5-2 keV corrected for the potential 
contribution from star formation; see §3.2. The X-ray data of HDF-oMD49 
are consistent with those expected for a heavily absorbed AGN, however, 
the low S/N X-ray spectrum alone provides limited diagnostics to distinguish 
between Compton-thin and Compton-thick absorption. We stress that since 
the data in this figure is normalised to the observed 4-8 keV flux it does not 
demonstrate the > 1 order of magnitude difference in flux between similarly 
powerful Compton-thin and Compton-thick AGNs; see Fig. 4. 



Krabbe etal.2001; Lutz etal.2004; Maiolino etal.2007; 
Horst etal. 2008). When combined with the observed X-ray 
emission, the rest-frame mid-IR luminosity can therefore pro- 
vide a diagnostic for the presence of heavy absorption. In this 
study we adopt the X-ray-6 pm luminosity relationship found 
by Lutz et al. (2004), which was determined by constrain- 
ing the AGN continuum component in large-aperture mid-IR 
spectra for a sample of nearby AGNs. We favour the Lutz 
et al. (2004) X-ray-mid-IR luminosity relationship over the 
small-aperture mid-IR imaging studies of Krabbe et al. (2001) 
and Horst et al. (2008) since HDF-oMD49 is unresolved by 
Spitzer. Furthermore, the Lutz et al. (2004) study covers a 
broader range in luminosity, and includes AGNs as luminous 
as those explored here. Indeed, we note that the X-ray-mid- 
IR luminosity ratio found for X-ray obscured but Compton- 
thin quasars (Sturm et al. 2006) is consistent with the Lutz 
et al. (2004) relationship, indicating that it is applicable for 
the objects explored in our study; see Fig. 4. However, the 
differences in the X-ray-mid-IR luminosity of these studies 
are small [the Lutz et al. (2004) luminosity ratio is within a 
factor of « 2-3 of that found by Krabbe etal. (2001) and 
Horst etal. (2008), assuming typical SEDs (i.e., NGC 1068; 
Mrk 231; M82) when converting to a common rest-frame 
wavelength of 6 pm] . 

Since it is now well established that the X-ray-to-optical lu- 
minosity ratio of optically selected AGNs is luminosity de- 
pendent (e.g., Vignali etal. 2003; Steffen etal. 2006), it is 
also interesting to explore whether the X-ray-mid-IR lumi- 
nosity ratio is luminosity dependent. Krabbe etal. (2001), 
Lutz etal. (2004), and Horst etal. (2008) did not find evi- 
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FIG. 4. — Rest-frame 2-10 keV luminosity versus 6 /jm AGN continuum luminosity. The observed and intrinsic properties of the z. ~ 2 AGNs selected from our 
study are shown; the X-ray data are converted from the observed-frame 0.5-2 keV band to the rest-frame 2-10 keV band assuming Y = 1.4. Five of the objects are 
identified using Spitzer-tRS and Chandra data (filled circle) and two are identified using Spitzer photometry and Chandra data (open circle); see Table 1. Objects 
indicated with a "G" lie in the GOODS-N field. The 6 fjm AGN luminosity is calculated from the best-fitting AGN component for the objects with Spitzer-tRS 
spectroscopy, and it is taken from the full SED for the objects with Spitzer photometry. The filled pentagon shows the intrinsic X-ray-6 /jm luminosity ratio 
of the X-ray obscured quasars studied by Sturm et al. (2006). The shaded region shows the range of intrinsic X-ray-mid-IR luminosity ratios found for AGNs 
in the local Universe (taken from Lutz et al. 2004), and the dotted lines indicate the observed rest-frame 2-10 ke\M) fim luminosity ratio expected for a typical 
AGN that is absorbed by Nn ~ 10 24 cm -2 (using the model presented in Alexander et al. 2005b); the dashed-dotted region shows the relationship that would be 
predicted assuming the luminosity-dependent Lx-Lmid-lR relationship of Maiolino et al. (2007). All of the plotted AGNs lie below the track for Nn « 10 24 cm -2 , 
suggesting that they are Compton thick. The dashed lines indicate the intrinsic 2-10 keV luminosities estimated from optical/near-IR spectroscopy (the individual 
emission-line luminosities are labelled) and provide good agreement with the intrinsic luminosities predicted from the mid-IR spectroscopy, confirming that all 
of the plotted AGNs are Compton thick. The average observed X-ray-IR luminosity ratio found for the mid-IR excess galaxies explored by Daddi et al. (2007a; 
open circle) is also shown, indicating that their average properties are consistent with those expected for Compton-thick AGNs about an order of magnitude less 
luminous than those explored here (assuming 100% of their objects are Compton-thick AGNs; see §4.4 for discussion); however, spectroscopic observations of 
individual objects (some of which are presented here) are required to provide robust confirmation of Compton-thick AGN activity. 



dence for a luminosity dependent X-ray-mid-IR luminosity 
ratio for local AGNs over w 3-5 orders of magnitude in lu- 
minosity. By comparison, the inference from the study of 
distant optically selected quasars by Maiolino et al. (2007) is 
that the X-ray-mid-IR luminosity ratio may be luminosity de- 
pendent. For example, solving Eqn. 1 and Eqn. 5 in Maiolino 
et al. (2007) gives an X-ray-mid-IR luminosity relationship of 
Lx x Lsjum s . However, there is considerable uncertainty in 
this result, and the derived X-ray-mid-IR luminosity relation- 
ship would be linear if a different X-ray-optical luminosity ra- 
tio from Steffen et al. (2006) is used (i.e., if Eqn. lb is adopted 
instead of Eqn. lc). In any case, the differences between the 
Lutz et al. (2004) and Maiolino et al. (2007) relationships are 
small (factor < 2-3) for the luminosities of the majority of 
the AGNs explored here (Lx < 3 x 10 44 erg s _1 ; see Fig. 4); 
however, we indicate the cases where our results would sig- 



nificantly change if the luminosity dependent version of the 
Maiolino et al. (2007) X-ray-mid-IR luminosity relationship 
is used. 

In Fig. 4 we show the rest-frame 2-10 keV luminosity ver- 
sus 6 fjm continuum luminosity of HDF-oMD49 and compare 
it to those found for AGN in the local Universe; we converted 
the rest-frame 6.7 /jm data used by Maiolino etal. (2007) to 
rest-frame 6 /jm assuming the SED of NGC 1068, which re- 
sults in small corrections (« 10%). The rest-frame 2-10 keV 
luminosity of HDF-oMD49 (L x = 1.4 x 10 42 erg s _1 ) is esti- 
mated from the observed-frame 0.5-2 keV luminosity (rest- 
frame 1.6-6.4 keV), for a spectral slope of V = 1.4. As- 
suming both the X-ray-mid-IR luminosity relationship found 
for local AGNs (Lutz et al. 2004), and the relationship de- 
rived from Maiolino et al. (2007), the intrinsic rest-frame 
2-10 keV luminosity of HDF-oMD49 is > 10 44 erg s" 1 , sug- 
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gesting that the X-ray absorption toward the AGN in HDF- 
oMD49 is Compton thick; see Fig. 4 and Table 1. However, 
to reduce further uncertainties on the X-ray-IR luminosity re- 
lationship due to, for example, the geometry and covering fac- 
tor of the gas and dust surrounding the X-ray emitting region, 
it is useful to have other measurements of the intrinsic AGN 
luminosity. 

The rest-frame UV spectrum of HDF-oMD49 shows strong 
AGN emission lines (CIV, CIH], Hell, Lyoc; see Fig. 2 of 
Steidel et al. 2002), which can be used to provide addi- 
tional estimates of the intrinsic AGN luminosity. We calcu- 
lated the emission-line fluxes using the Steidel et al. (2002) 
optical spectrum and integrating the flux over the emission- 
line profiles; see Reddy et al. (2006) for more details. 
Adopting the average emission-line ratios of Netzer et al. 
(2006), to calculate the expected [OIII]A,5007 luminosity, 
and the [OIII]X5007-X-ray flux ratio found by Mulchaey 
et al. (1994), the range of predicted intrinsic rest-frame 
2-10 keV luminosity from the different emission lines is 
(1.8^1.8) xlO 44 erg s -1 . We get similar intrinsic X-ray lu- 
minosities if we use the [OIII]A,5007-X-ray flux ratio found 
by Alonso-Herrero etal. (1997) and we get higher intrinsic 
X-ray luminosities (up to an order of magnitude higher for 
the most luminous AGNs) if we use the Netzer et al. (2006) 
[OIII]A,5007-X-ray flux ratio. We have not corrected these 
emission-line fluxes for extinction or contamination from star 
formation, which are poorly constrained at rest-frame UV 
wavelengths. However, we note that the Mulchaey et al. 
(1994), Alonso-Herrero etal. (1997), and Netzer etal. (2006) 
correlations were also derived using uncorrected emission- 
line fluxes. The uncertainties in both the emission-line ratios 
and the [OIII]A,5007-X-ray flux ratio are « 2 (see Table 1 of 
Netzer etal. 2006 and Table 2 of Mulchaey etal. 1994), giv- 
ing a combined uncertainty of a factor of « 3, comparable to 
the uncertainty in the intrinsic X-ray luminosity based on the 
rest-frame 6 luminosity. The intrinsic rest-fame 2-10 keV 
luminosities predicted from both the optical and the mid-IR 
spectroscopy are in good agreement, providing compelling 
evidence for Compton-thick absorption of an intrinsically lu- 
minous quasar in HDF-oMD49 (t-2-10 keV ~ 3 x 10 44 erg s -1 ; 
N H > 10 24 cm" 2 ). 

An alternative scenario for the weak X-ray emission from 
HDF-oMD49 is that the AGN has recently faded (e.g., 
Guainazzi etal. 1998). However, since the light-crossing time 
of the mid-IR emitting region is likely to be short (i.e., on 
the basis of the dust-sublimation radius it will be w 1-2 yrs; 
e.g., see Table 1 of Granato, Danese, & Franceschini 1997 and 
Fig. 30 of Suganuma et al. 2006), this scenario is only tenable 
for HDF-oMD49 if we have caught the AGN during a very 
brief transitionary stage where the X-ray emission has faded 
but the mid-IR emission is still strong. The good agreement 
between the fluxes measured from the Spitzer-MlPS photom- 
etry and the Spitzer-IRS spectroscopy, taken 2 years apart, 
suggests that this scenario is unlikely. 

4. DISCUSSION 

We have used Spitzer-IRS spectroscopy, 3.6-70 pm pho- 
tometry, optical spectroscopy, and ultra-deep 2 Ms Chandra 
data to show that HDF-oMD49 hosts an intrinsically lumi- 
nous AGN that is heavily obscured by Compton-thick mate- 
rial (L 2 -io keV ~ 3 x 10 44 erg s _1 ; N H > 10 24 cm~ 2 ); i.e., it is 
a Compton-thick quasar. The S/N of the X-ray data alone was 
too poor to show conclusively that the AGN is obscured by 
Compton-thick material. However, when combined with rest- 



frame 6 pm and optical emission-line luminosity constraints, 
the evidence for Compton-thick absorption in this X-ray faint 
source is compelling. Additional support for this interpre- 
tation comes from the similarity of the mid-IR spectrum of 
HDF-oMD49 and the composite mid-IR spectrum of X-ray 
obscured but Compton-thin quasars; see Fig. lb. 

On the basis of the comparatively weak ultra-hard 
4-8 keV emission from HDF-oMD49, even the rest-frame 
12.8-25.7 keV emission appears to be heavily suppressed. 
For example, assuming an X-ray spectral slope of T = 2.0, 
the predicted 12.8-25.7 keV luminosity based on the intrinsic 
rest-frame 2-10 keV luminosity is a factor of « 20 higher than 
the observed 12.8-25.7 keV luminosity (m 1.3 x 10 44 erg s _1 
versus w 6 x 10 42 erg s _1 ; see §2.2). It is difficult to relate 
this factor of 20 suppression to an accurate estimate of the 
absorbing column density since the observed X-ray emission 
in the 12.8-25.7 keV band could be dominated by reflected 
and scattered components (e.g., NGC 1068; Matt etal. 1997; 
see Fig. 8.1 in Comastri 2004). However, on the basis of 
the PLCABS model (Yaqoob 1997) in XSPEC, which considers 
the effects of electron scattering and photoelectric absorption 
through an absorbing medium, this amount of suppression in 
the 12.8-25.7 keV band would correspond to an absorbing 
column density of Nn ~ 4 x 10 24 cm . Due to the absence 
of a reflection component in this model, and the simplify- 
ing assumption of a spherical geometry, this absorbing col- 
umn density should only be considered representative (e.g., 
see Fig. 8.1 in Comastri 2004 for an alternative model). 

In this final section we adopt the approach that we used 
in §3 to identify other distant Compton-thick AGNs using 
data obtained from the literature. Utilizing these constraints 
we then directly estimate the space density of Compton-thick 
quasars at z ~ 2 and compare our results to current observa- 
tional and theoretical constraints. 

4. 1 . The Identification of other Distant Compton-thick AGNs 

From a search of the published Spitzer-IRS spectroscopic 
data, there are four good candidate z ~ 2 Compton-thick 
AGNs that have mid-IR spectroscopy and the essential sen- 
sitive X-ray constraints required to identify Compton-thick 
AGN activity: FSC 10214+4714, SW J104406+583954, 
SW J104409+585224, and SMM J123600+621047. 

FSC 10214+4714 is a strongly lensed Seyfert 2 galaxy at 
Z = 2.285 which produces luminous 6 pm AGN continuum 
emission (Teplitz et al. 2006) but is only weakly detected at 
X-ray energies, despite an effective lensing-corrected Chan- 
dra exposure of m 2 Ms (Alexander et al. 2005a). In Fig. 4 
we plot FSC 10214 and provide secondary support for the 
intrinsic luminosity of the AGN using the [OIII]A.5007 lumi- 
nosity from Alexander et al. (2005a) and the [OIH]A,5007- 
X-ray flux ratio of Mulchaey et al. (1994); see Table 1. 
The combination of the X-ray data and optical-mid-IR spec- 
troscopy confirms that FSC 10214 hosts a Compton-thick 
quasar (L2-i0keV ~ 2 x 10 44 erg s -1 ), even given the consid- 
erable scatter in the relationships. 

SW J104406+583954 and SW J104409+585224 are X-ray 
weak candidate Compton-thick AGNs identified in the Chan- 
dra/S WIRE survey (Polletta et al. 2006) that lie at at z = 2.430 
and z = 2.540, respectively. The optical and mid-IR spectra 
of both sources show clear signatures of AGN activity (Pol- 
letta etal. 2006, 2008; Weedman etal. 2006), which we use 
to determine the intrinsic luminosities of the central sources; 
see Table 1. In Fig. 4 we plot their X-ray-6 pm luminosi- 
ties, showing that the properties of both objects are consis- 
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tent with those expected for a luminous Compton-thick quasar 
(^2-iOkeV > 10 45 erg s _1 ); however, we note that the evidence 
is weaker for SW 104406 if the luminosity dependent X- 
ray-to-mid-IR luminosity relationship inferred from Maiolino 
et al. (2007) is used. 

SMM J 123600+62 1047 is a z = 2.002 submillimeter emit- 
ting galaxy (SMG) with a mid-IR-bright AGN (Pope et al. 
2008) that is undetected in the 2 Ms CDF-N data (Alexan- 
der et al. 2003a); we ran WAVDETECT using a false-positive 
probability threshold of 10~ 5 but did not detect X-ray emis- 
sion from this source (see Table 1 of Alexander etal. 2005b). 
The X-ray-6 pm luminosity ratio indicates that it hosts a 
Compton-thick AGN (see Fig. 4), although optical and near- 
IR spectroscopic observations do not reveal the signatures of 
AGN activity (e.g., Swinbank et al. 2004; Chapman et al. 
2005), leaving some uncertainty on the intrinsic luminosity of 
the AGN. However, since the mid-IR luminosity of the AGN 
in SMMJ 123600 is about an order of magnitude larger than 
that found in typical SMGs (e.g., Lutz et al. 2005; Menendez- 
Delmestre et al. 2007; Valiante et al. 2007; Pope et al. 
2008), the w 10 times larger intrinsic X-ray luminosity de- 
rived here seems plausible (Lx.6^m ~ 5 x 10 44 erg s , i.e., 
compare to Alexander et al. 2005b); see Table 1. In Alexan- 
der et al. (2005c), it was argued that SMGs represent a rapid 
black-hole growth phase before the onset of optically lumi- 
nous quasar activity. The discovery of this one mid-IR lu- 
minous X-ray undetected SMG in the CDF-N field increases 
current constraints on the integrated black-hole growth den- 
sity in SMGs by a factor of w 2 (see Fig. 2 of Alexander et al. 
2005c). 

Other studies have used optical spectroscopy to identify 
AGNs that are mid-IR bright and X-ray weak/undetected but 
lack Spitzer-IRS spectroscopy. Reddy et al. (2006) found two 
z f» 2.4-2.5 optically identified AGNs (BX160; BX1637) that 
are X-ray undetected in the 2 Ms CDF-N survey and have 
power-law like IR SEDs [the 24 pm fluxes of these sources 
(fiApm ~ 100-140 /jJy) are too faint for good S/N Spitzer- 
IRS spectroscopy]; see Table 1. Following our approach for 
HDF-oMD49 (see §2.2), we ran WAVDETECT using a false- 
positive probability threshold of 10~ 5 but did not detect X-ray 
emission from either source; 3 O X-ray upper limits are cal- 
culated following §3.4.1 in Alexander et al. (2003a). We plot 
the X-ray-to-6 pm luminosity ratios of these objects in Fig. 4. 
The mid-IR data, when combined with the optical emission- 
line luminosities, indicates that these sources are likely to 
be Compton-thick quasars; we determined their emission-line 
properties using the same procedure as for HDF-oMD49 (see 
§3.3). The evidence would be weaker for BX1637 if the lu- 
minosity dependent X-ray-mid-IR luminosity ratio inferrred 
from Maiolino et al. (2007) is used, although we note that the 
X-ray-6 pm luminosity ratio for this source is an upper limit. 

We also note that other studies have selected candi- 
date Compton-thick AGNs using a combination of multi- 
wavelength data with sensitive X-ray constraints (e.g., Donley 
et al. 2005, 2007; Alonso-Herrero 2006; Daddi et al. 2007a; 
Martiriez-Sansigre et al. 2007; Fiore et al. 2008); many other 
studies have identified obscured AGNs using mid-IR spec- 
troscopy but lack the essential X-ray data to provide a case for 
Compton-thick absorption. However, none of these studies 
has optical or mid-IR spectroscopy to provide accurate mea- 
surements of the intrinsic AGN luminosities. 
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FIG. 5. — Rest-frame 2-10 keV luminosity versus spectroscopic redshift 
for the X-ray detected sources in the GOODS-N region of the CDF-N sur- 
vey and the seven Compton-thick AGNs identified here; the four Compton- 
thick AGNs that lie in the GOODS-N field are indicated with a "G". The 
dotted curve shows the aim-point sensitivity limit for the CDF-N survey; see 
Fig. 18 of Alexander et al. (2003a) for the drop in Chandra sensitivity for 
off-axis sources. The source redshifts are predominantly taken from Barger 
et al. (2003). The dashed lines connect the observed X-ray luminosity of 
the Compton-thick AGNs to their intrinsic luminosities; the intrinsic lumi- 
nosities correspond to the mean luminosity estimated from the combination 
of the 6 fim luminosity and the average emission-line luminosity given in 
Table 1. Although the four Compton-thick AGNs that lie in the GOODS-N 
region are either X-ray weak or undetected, they are amongst the most intrin- 
sically luminous AGNs in the GOODS-N field. 



4.2. The Properties of Distant Compton-thick AGNs 

The intrinsic X-ray-mid-IR continuum ratios, inferred us- 
ing the rest-frame UV emission-line luminosities of the six 
Compton-thick quasars with optical AGN signatures, are con- 
sistent with those found for AGNs in the local Universe; 
see Fig. 4. This is in agreement with the results of Sturm 
et al. (2006), who found that the average intrinsic X-ray-6 pm 
luminosity of distant X-ray obscured quasars are similar to 
those found for local AGNs, suggesting that the objects ex- 
plored here are the extreme Compton-thick end of the X-ray 
obscured quasar population. We also showed in §3.1 and 
Fig. lb that the Spitzer-IRS spectrum of HDF-oMD49 is sim- 
ilar to the composite mid-IR spectrum produced by Sturm 
et al. (2006), providing a particularly strong case for a quasar 
that is absorbed in the X-ray band by Compton-thick material 
in this source. These pieces of evidence also suggest that the 
accretion properties of distant optically identifiable AGNs are 
similar to those observed in local AGNs, in agreement with 
previous studies that have explored the X-ray-optical and X- 
ray-infrared properties of optically selected AGNs (e.g., Stef- 
fen et al. 2006; Jiang et al. 2006). 

In Fig. 5 we show the observed and intrinsic X-ray lumi- 
nosities versus redshift of the seven Compton-thick quasars 
identified here. Four of the Compton-thick quasars lie in the 
GOODS-N region (w 160 arcrrnV) of the CDF-N and, despite 
being X-ray weak/undetected, they are amongst the most in- 
trinsically luminous AGNs identified in this field. This shows 
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FIG. 6. — Selection of distant Compton-thick AGNs. (a) ratio of mid-IR-based to UV-based star-formation rate (SFR) versus 8 /im luminosity for the z ~ 2 
galaxies studied by Daddi etal. (2007 a,b) and the spectroscopically identified Compton-thick quasars; the open and filled symbols correspond to objects with 
photometric and spectroscopic redshifts, respectively. This figure has been adapted from Fig 2 of Daddi et al. (2007a). (b) mid-IR-optical flux ratio versus R — K 
color for the galaxies investigated by Fiore et al. (2008) and Daddi et al. (2()()7a,b) and the spectroscopically identified Compton-thick quasars; see also Fig. 3 of 
Fiore etal. (2008). FSC 10214+4724 is not plotted on these figures since there is some uncertainty on the lensing magnification of the star-forming component. 
The dashed lines indicate the selection criteria for the Daddi etal. (2007a) and Fiore etal. (2008) studies. Our spectroscopically identified Compton-thick quasars 
are extreme examples of the mid-IR excess galaxy population of Daddi et al. (2007a) but none of them formally match the Fiore et al. (2008) selection criteria, 
which is probably due to the requirement for an optical signature of AGN activity in our sample. This implies that there could be many more distant Compton-thick 
AGNs than that suggested by either study alone. 



that even in the deepest X-ray survey currently available, we 
are only able to directly probe the peak of the Compton- 
thick AGN population at z ~ 2. The two luminous Compton- 
thick quasars were identified in the wider-area SWIRE survey 
(w 0.6 deg 2 ) and are amongst the most luminous obscured 
AGNs known. 

Four of the seven Compton-thick quasars are detected in the 
X-ray band, most typically at > 2 keV. Three of these objects 
have hard X-ray spectra (T < 1; HDF-oMD49; SW 104406; 
SW 104409; Polletta etal. 2006) and one object has a com- 
paratively softX-ray spectrum (T« 1.6; FSC 10214; Alexan- 
der et al. 2005c). The comparatively soft X-ray spectrum 
of FSC 10214 might be dominated by scattered emission 
and star formation from the host galaxy, as often found in 
Compton-thick AGNs (e.g., Bassani et al. 1999; Comastri 
2004; Guainazzi et al. 2005). We can also place basic X-ray 
spectral constraints on the three X-ray undetected Compton- 
thick quasars by stacking their X-ray data; all of these objects 
lie in the GOODS-N field where the X-ray data are partic- 
ularly sensitive. Using the X-ray stacking code adopted by 
Worsley etal. (2005), we obtain » 3 a detections in the 0.5- 
8 keV and 2-8 keV bands, and have a«2o constraint in 
the 0.5-2 keV band; the significance of these stacking results 
were calculated using 10,000 Monte Carlo trials. The stacked 
data corresponds to fluxes (and a 3 a upper limit) of 2.8 x 
10~ 16 erg s- 1 cm- 2 (0.5-8 keV), < 2.3 x 10~ 17 erg s" 1 cm" 2 
(0.5-2 keV), and 2.7 x 10~ 16 erg s" 1 cm~ 2 (2-8 keV). The 
average rest-frame 2-10 keV luminosity constraint for these 
three objects derived from the observed-frame 0.5-2 keV 
flux is < 10 42 erg s -1 . The derived X-ray spectral slope of 



r < 0.3 unambiguously confirms that heavily obscured AGNs 
are present in these sources (e.g., see Fig. 2a in Alexander 
etal. 2005b). 

4.3. The Selection of Distant Compton-thick AGNs 

The seven Compton-thick quasars investigated here were 
selected on the basis of unambiguous evidence for AGN ac- 
tivity from optical and/or mid-IR spectropscopy. The advan- 
tage with this approach is that we have been able to esti- 
mate the intrinsic AGN properties of these objects using the 
optical-mid-IR spectroscopy. However, the requirement for 
spectroscopic evidence of AGN activity has led to a compar- 
atively small sample of objects. Nevertheless, we can provide 
some insight into the completeness of Compton-thick AGN 
selection by comparing the properties of our reliably identi- 
fied Compton-thick quasars to the larger samples of candidate 
Compton-thick AGNs that have been photometrically identi- 
fied in other studies. Arguably, the most complete identifica- 
tion studies of distant Compton-thick AGNs to date are Daddi 
etal. (2007a,b) and Fiore et al. (2008). 

Daddi et al. (2007a,b) identified a large population of X-ray 
undetected z ~ 2 galaxies with an excess of mid-IR emission 
over that expected from star formation (as predicted using 
the dust-extinction corrected UV luminosity). X-ray stack- 
ing analyses revealed a hard X-ray spectral slope from these 
objects (r = 0.8), comparable to that of the Compton-thick 
quasars and unambiguously identifying the presence of heav- 
ily obscured AGNs; see §4.2. In Fig. 6a we plot the ratio 
of mid-IR based to UV-based star-formation rate (SFR) ver- 
sus 8 /jm luminosity of the Compton-thick quasars and com- 
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pare them to Daddi et al. (2007a,b); we have not calculated the 
properties of FSC 10214 due to the significant uncertainties 
in the lensing magnification of the star-forming regions. The 
Compton-thick quasars lie significantly above the mid-IR ex- 
cess threshold defined by Daddi et al. (2007a), indicating that 
they are extreme examples of the mid-IR excess galaxy pop- 
ulation. This analysis shows reliably that at least a subsample 
of the mid-IR excess galaxy population host Compton-thick 
AGNs. 

Fiore et al. (2008) took a complementary approach to that of 
Daddi et al. (2007 a,b) and selected X-ray undetected objects 
with extreme mid-IR-to-optical flux ratios (/24 A ,m/ '/r > 1000) 
and red optical colors (R—K> 4.5). X-ray stacking analyses 
of these objects revealed a hard X-ray spectrum with a slope 
similar to that found for the mid-IR excess galaxies (r« 1 .0, 
calculated from their stacked X-ray count rates). In Fig. 6b 
we plot the mid-IR-to-optical flux ratio versus R — K color of 
the Compton-thick quasars and compare them to the candi- 
date Compton-thick AGNs studied by Fiore et al. (2008) and 
Daddi et al. (2007a); as before, we did not calculate the prop- 
erties of FSC 10214. Formally, none of the Compton-thick 
quasars match the candidate Compton-thick AGN selection 
criteria of Fiore et al. (2008), although two of the objects lie 
close. Furthermore, only w 10% of the mid-IR excess galax- 
ies of Daddi etal. (2007a) match the Fiore etal.(2008) cri- 
teria. These analyses show that the candidate Compton-thick 
AGN criteria of Daddi et al. (2007a) and Fiore et al. (2008) are 
mutually exclusive, suggesting that the distant X-ray unde- 
tected AGN population may be larger than either study has 
suggested. The requirement in our study for optical spectro- 
scopic redshifts and the optical identification of AGN activity 
(with the exception of SMM J 123600) probably causes a bias 
against selecting objects with red R —K colors, which by def- 
inition will be optically faint. 

4.4. The Ubiquity of Distant Compton-thick AGNs 

Our seven Compton-thick quasars do not comprise a com- 
plete sample. However, four of the objects lie in the 
GOODS-N field, providing basic constraints on the space 
density of Compton-thick quasars at z ~ 2-2.5. On the ba- 
sis of the comoving volume at z = 2-2.5 in a region the size 
of the GOODS-N field, we estimate a comoving space density 
for Compton-thick quasars withL2-i0keV ~ 1 44 — 1 45 ergs -1 
of <t> w (0.7-2.5) x 10~ 5 Mpc~ 3 , where the range only reflects 
the uncertainty due to small-number statistics (e.g., Gehrels 
1986). 

In Fig. 7 we plot the Compton-thick quasar space density 
and compare it to that found in other studies. Our results 
suggest that the space density of Compton-thick quasars is 
w 1-5 times higher than that of comparably luminous un- 
obscured quasars at z ~ 2-2.5 ($«5x 10~ 6 Mpc~ 3 ; e.g., 
Hasinger et al. 2005). Since the Gilli et al. (2007) model of the 
X-ray background constrains the number of Compton-thick 
quasars to be the same as the number of unobscured quasars, 
our Compton-thick quasar space density is also m 1-5 times 
higher than the Compton-thick quasar predictions of Gilli 
etal. (2007). As we mention in §4.1 and show in Fig. 4, the 
evidence for a Compton-thick quasar in BX1637 is weaker 
if the X-ray-mid-IR luminosity ratio inferred by Maiolino 
et al. (2007) is used. However, we note that a significantly 
more important issue is likely to be the effect of cosmic vari- 
ance since we have only identified Compton-thick quasars in 
a small region over a narrow redshift range. For example, 
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FIG. 7. — Space density of Compton-thick AGNs. The filled square 
corresponds to the results determined here for spectroscopically identified 
Compton-thick quasars with L2-i0keV ~ 10 44 -10 45 erg s _I that lie in the 
GOODS-N region of the CDF-N field. The open squares correspond to the 
results obtained by Daddi et al. (2007a) and Fiore et al. (2008) using photo- 
metrically identified objects and X-ray stacking analyses; the plotted values 
assume that 100% of their objects are Compton-thick AGNs and they are 
therefore upper limits (see §4.4 for discussion). The dashed lines are the 
predicted space densities of Compton-thick AGNs for a range of X-ray lumi- 
nosity lower limits from the Gilli et al. (2007) XRB model. The space density 
of the Compton-thick quasars robustly identified here is a factor 1-5 higher 
than that found for comparably luminous unobscured quasars, indicating that 
Compton-thick accretion is ubiquitious at z ~ 2; see §4.4 for a more detailed 
discussion. 



within the GOODS-N field, there are also four X-ray unob- 
scured quasars that lie at z ~ 2-2.5, indicating that the true 
Compton-thick-unobscured quasar ratio may be toward the 
lower end of the value given above. Nevertheless, our results 
clearly indicate that a large fraction of the growth of black 
holes was obscured by Compton-thick material, in general 
agreement with results based on less-reliable photometrically 
identified objects (e.g., Daddi etal. 2007a; Martmez-Sansigre 
et al. 2007; Fiore etal. 2008), the predictions made by theo- 
retical models (e.g., Marconi et al. 2004; Treister et al. 2006), 
and the Compton-thick AGN fraction inferred by the study of 
Maiolino et al. (2007; Fig. 7a). 

We can compare our estimated space density to those deter- 
mined from the photometric selection of candidate Compton- 
thick AGNs by Daddi et al. (2007a) and Fiore et al. (2008). 
On the basis of the X-ray data and the X-ray-6 /jm lu- 
minosity, Daddi et al. (2007a) argued that their candidate 
Compton-thick AGNs have an average intrinsic luminosity 
of L2-i0keV ~ (1-4) xlO 43 erg s . The intrinsic lumi- 
nosities of these mid-IR excess galaxies are about an order 
of magnitude lower than the Compton-thick quasars iden- 
tified here, and consequently their space density is higher 
(4> w 2.6 x 10~ 4 Mpc~ 3 ; see Fig. 7). Fiore etal. (2008) es- 
timated the intrinsic X-ray luminosities of their X-ray unde- 
tected candidate Compton-thick AGNs using the 5.8 lu- 
minosities, finding L2-i0keV > 10 43 erg s . Taking into ac- 
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count the estimated redshift range of their objects (z ~ 1.2- 
2.6), we calculate a space density of <t> w 7 x 10~ 5 Mpc~ 3 ; as 
mentioned in Fiore et al. (2008), this space density is consis- 
tent with that predicted for Compton-thick AGNs using the 
XRB model from Gilli etal. (2007). Taken at face value, 
these results suggest that lower-luminosity Compton-thick 
AGNs were « 5-20 times more common than Compton-thick 
quasars at z « 2. However, since the primary evidence for 
Compton-thick AGN activity in these studies was from X- 
ray stacking analyses, there could be significant numbers 
of non-Compton-thick AGNs in their samples (i.e., either 
Compton-thin AGNs or star-forming galaxies). Both Daddi 
et al. (2007a) and Fiore et al. (2008) asssessed the "contamina- 
tion" from star-forming galaxies to their stacked X-ray results 
and concluded that > 50% and « 80 ± 15% of the objects 
are heavily obscured AGNs, respectively. However, these per- 
centages refer to the fraction of photometrically identified ob- 
jects that are heavily obscured AGNs and direct optical-mid- 
IR spectroscopic identification will be required to determine 
what fraction of these AGNs are obscured by Compton-thick 
material. 

We can also compare our results to those determined 
for more luminous Compton-thick quasars. Polletta 
etal. (2006) identified five candidate Compton-thick quasars 
at z~ 1.4-2.6 with L2-i0keV > 10 45 erg s _1 on the basis 
of the Chandra data in the SWIRE survey, giving a space 
density of <t> « 5.7 x 10~ 7 Mpc~ 3 . Two of these objects 
have been presented here as robust Compton-thick quasars 
(SW J104406 and SW J104409; see Table 1). However, the 
other three candidate Compton-thick quasars are all unde- 
tected at 24 pm (jy^im < 230 piy; see Table 2 of Polletta 
et al. 2006). This suggests that either the intrinsic AGN lu- 
minosities of these Compton-thick quasars are more modest 
(^2-iOkeV < 10 45 erg s~ ! ) than those estimated by Polletta 
et al. (2006) or that the X-ray-mid-IR relationship shown in 
Fig. 4 does not hold for these luminous objects. Martmez- 
Sansigre et al.(2007) identified ten Compton-thick quasars 
at z ~ 1.7-4.2 with an average bolometric luminosity of 
Lbol ~ 10 47 erg s -1 , giving a space density of 4> w 3.9 x 
10~ 8 Mpc~ 3 . On the basis of the bolometric conversion 
of Elvis etal. (1994), which Martmez-Sansigre et al.(2007) 
also adopted, the estimated average X-ray luminosities are 
^2-iOkeV ~ 3 x 10 45 erg s -1 for these sources, suggesting 
that they are very luminous AGNs. However, since two of 
their objects are Compton-thin quasars that are detected at X- 
ray energies, we can directly estimate the intrinsic 2-10 keV 
luminosities of the candidate Compton-thick quasars using 
the average X-ray-6 luminosity ratio of the Compton- 
thin quasars. Taking this approach and using the median 
24 /jm flux for the sample, the typical intrinsic X-ray lumi- 
nosities of these quasars would be more modest (L2-10 keV ~ 
3 x 10 44 erg s _1 ), casting significant doubt on whether they 
are obscured by Compton-thick material; the rest-frame 6 pm 
fluxes are derived from the observed 24 pm fluxes using 
the NGC 1068 SED to apply small -corrections. Direct 
spectroscopic identification of these photometrically classi- 
fied Compton-thick quasars is required to place more direct 
constraints, which should be possible with current instrumen- 
tation. 

Lastly, we compare our results to the Compton-thick AGNs 
identified by Tozzi etal. (2006) using X-ray spectral anal- 
ysis of X-ray detected AGNs in the Chandra Deep Field- 
South survey. Fourteen of the objects investigated by Tozzi 



et al. (2006) had X-ray spectra that were better fitted by a 
pure reflection model than an absorbed power-law model, and 
were consequently classified as Compton-thick AGNs; de- 
rived space densities are 4> w 3 x 10~ 5 Mpc~ 3 at z « 1 and 
4> w 8 x 10~ 6 Mpc~ 3 at z ~ 2. While the Tozzi et al. (2006) 
study provides the most quantitative X-ray identification of 
distant Compton-thick AGNs to date, it is also possible for 
Compton-thin AGNs to have a pure reflection component (i.e., 
if the power-law component has temporarily decreased), and 
other pieces of evidence are required to confirm that these are 
Compton-thick AGNs. Two of the seven objects with spectro- 
scopic redshifts have strong Fe Ka emission lines, confirming 
that they are likely to be Compton thick, but supporting evi- 
dence for Compton-thick absorption is lacking in the other 
objects. 

4.5. Prospects for Improved Observational Constraints 

The current study provides diagnostics to identify distant 
Compton-thick AGNs but is limited in source statistics, cov- 
ers a narrow redshift range, and is restricted to the bright- 
est AGNs at z « 2. Future studies could focus on obtain- 
ing optical-mid-IR spectroscopy of lower-luminosity z ~ 2 
candidate Compton-thick AGNs, as well as identifying ob- 
jects at lower and higher redshifts. Assuming the NGC 1068 
template used here and the average X-ray-mid-IR relation- 
ship shown in Fig. 4, an AGN with L2-I0kev ~ 10 43 erg s _1 
will have 24 pm fluxes of w 2 pjy, w 15 piy, « 0.3 mJy, 
and ~ 5 mJy at z « 4, z « 2, z w 0.7, and z w 0.2, respec- 
tively; these fluxes will typically scale linearly with lumi- 
nosity (see §3.3). Mid-IR bright objects (/24/jm > 0.2 mJy) 
should be identifiable with Spitzer-IRS but mid-IR spectro- 
scopic identification of fainter objects will need to wait un- 
til the launch of the James Webb Space Telescope. Optical 
and near-IR spectroscopy may be easier to obtain for many 
of the candidate Compton-thick AGNs, although only if the 
AGNs have identifiable emission lines, which could be weak 
in low-luminosity and dust-reddened AGNs (e.g., Caccianiga 
etal. 2007; see §4.3). 

Ultimately, deep X-ray data are required to search directly 
for the presence of Compton-thick absorption in individual 
objects using high S/N X-ray spectroscopy. For the major- 
ity of the objects investigated here, this will require ^> 2 Ms 
exposures with current X-ray telescopes; see Fig. 5. How- 
ever, due to its large lensing boost it is possible to provide 
reasonable X-ray spectral constraints for FSC 10214+4724 in 
a modest Chandra exposure (e.g., a 100 ks Chandra exposure 
would yield almost 100 X-ray counts, > 10 times more than 
found for HDF-oMD49; see Fig. 3). Large improvements 
in the X-ray spectral constraints for distant Compton-thick 
AGNs will need to wait for the next generation of X-ray ob- 
servatories (e.g., Constellation-X, NuSTAR, Simbol-X,XEUS), 
where the foci are high-resolution spectroscopy and/or high- 
energy imaging (> 10 KeV) of faint X-ray sources. For ex- 
ample, a 1 Ms XEUS exposure will yield w 5000 counts for 
HDF-oMD49, allowing detailed X-ray spectral analyses and 
providing the potential to search for changes in the X-ray 
emission due to column density variations, giving direct in- 
sight into the environment around the black hole of a distant 
Compton-thick AGN (e.g., Risaliti et al. 2002). 

5. CONCLUSIONS 

We have presented Spitzer-IRS spectroscopy and 
3.6-70 Ami photometry of HDF-oMD49, a z = 2.211 
optically identified AGN that is formally undetected in 
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the 2 Ms CDF-N observation. From a combination of 
optical-mid-IR spectroscopy and X-ray data, we have 
shown that HDF-oMD49 hosts an intrinsically lumi- 
nous quasar that is obscured by Compton-thick material 
(£2-10 kev « 3 x 10 44 erg s" 1 ; N H > 10 24 cm" 2 ). We 
selected six further z ~ 2 AGNs from the literature (four 
with Spitzer-IRS spectroscopy) and used the same X-ray- 
optical-mid-IR diagnostics applied to HDF-oMD49 to show 
that these objects are also likely to be Compton-thick quasars 
with L2-i0keV > 10 44 erg s _1 . We demonstrated that these 
Compton-thick quasars would be classified as mid-IR excess 
galaxies, on the basis of the Daddi et al. (2007a) definition, 
providing the first spectroscopic confirmation of Compton- 
thick AGN activity in a subsample of these z ~ 2 mid-IR 
bright galaxies. Four of these Compton-thick quasars he in 
the GOODS-N field, and we used these objects to constrain 
the space density of distant Compton-thick quasars, finding 
4> w (0.7-2.5) xl0~ 5 Mpc- 3 at z « 2-2.5. On the basis 
of our results, Compton-thick quasars were as ubiquitious 
as unobscured quasars at z « 2-2.5, implying that a large 
fraction of the growth of supermassive black holes must have 
been obscured by Compton-thick material. 
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